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ABSTRACT 

This work presents a new approach to gyro calibration where, in addition to being used for 
computing attitude that is needed in the calibration process, the gyro outputs are also used as measurements 
in a Kalman filter. This work also presents an algorithm for calibrating a quadruplet rather than the 
customary triad gyro set. In particular, a new misalignment error model is derived for this case. The new 
calibration algorithm is applied to the EOS-AQUA satellite gyros. The effectiveness of the new algorithm is 
demonstrated through simulations. 


INTRODUCTION 

Gyro calibration as well as calibration of other instruments includes two stages. In the first stage 
the instrument error parameters are estimated. During the second stage those errors are continuously 
removed from the gyro readings. In the classical approach to gyro calibration, the gyro outputs are used to 
maintain or compute body orientation rather than being used as measurements in the context of filtering. In 
inertial navigation, for example (ref. 1), gyro errors cause erroneous computation of velocity and position, 
and then when the latter are compared to measured velocity and position, a great portion of the computed 
velocity and position errors can be determined. The latter errors are then fed into a Kalman filter (KF) that 
uses the Inertial Navigation System (INS) error model to infer on the gyro errors. Similarly, when applying 
the classical approach to spacecraft (SC) attitude determination, the gyro outputs are used to compute the 
attitude and then attitude measurements (refs. 2, 3) are used to determine the attitude errors, which again 
using a KF, indicates what the gyro errors are. 

In the approach adopted in this work, the gyro outputs are used as angular rate measurements and 
are compared to estimated angular rate measurements. However, this approach requires the knowledge of 
the angular rate. In the past (ref. 4), the estimated angular rate was computed in a rather simplistic way 
assuming basically that the rate was constant. In the present work, the estimated angular rate is derived 
using a KF whose input can be any kind of attitude measurement; therefore, the angular rate experienced by 
the SC can be continuously changing, and yet a good estimate of the rate, necessary for calibration, can be 
obtained. 


The calibration algorithm presented in this work was derived for a set of quadruplet gyros. This 
required the derivation of a new error model, particularly for the gyro misalignments. The new calibration 
algorithm was applied to the gyro package of the EOS-AQUA satellite. The latter consists of four gyros, 
which are given the task of measuring the three components of the SC angular velocity vector resolved in 
the body Cartesian coordinates. 

In the next section the gyro error model is derived. The section that follows presents an algorithm 
for computing the calibration parameters when the rate is known, and then in the section that follows we 
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present the same when the rates are unknown. In the following section we present the compensation 
procedure that needs to take place to complete the calibration process, and in the subsequent section we 
present simulation results. Finally, in the last section, the conclusions are presented. 


GYRO ERROR MODEL 

The gyro errors that are considered in this work are: misalignment, scale factor error, and bias 
(constant drift rate). The gyro error model is basically a linear model, which associates small error sources 
to the gyro outputs. Due to the linearity of the model we can compute the contribution of each error source 
independently, and then sum up all the contributions into one linear model. 

We start the description of the error model, by deriving the expression for the gyro misalignments. 


Actual 
Direction 
Gyro No. 


Assumed 
(nominal) 
Direction of 
Gyro No. j 


Y, co. 


X, co 


X 



Fig. 1 : The geometry of the Assumed and the Actual Direction of the Gyro Input Axis 


Misalignment Model 

The assumed direction of the sensitive axis of gyro j, which is one of the four gyros, is presented 
in Fig. 1 where the body coordinate axes are also presented and are denoted by X, Y, and Z. The orientation 
of this gyro is expressed by a vector of unit length in the direction of the gyro sensitive axis. The direction 
of this unit vector in the body coordinates is expressed by its three direction cosines, which are identical to 
its components when the unit vector is resolved in the body coordinates. These components are 
Cjj, Cj 2 , and C j3 . Being direction cosines, or equivalently, components of a unit vector, the sum of their 
squares adds up to 1 ; that is, 

C jl + C j2 + C j3 = 1 C 1 ) 

The rate that this gyro reads is the projection of the angular velocity vector on this unit vector. If 
we express the angular velocity vector in the body coordinates, where its components are co x , © y , and (D z , 
then this projection is given by 

lj ■ w = Cj,a) x + c j2 w y + c j3 © z (2.a) 
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where lj is the unit vector along the j— gyro sensitive axis, and CO is the angular rate vector. The 
nominal (error-less) reading of this gyro is then 

r» x ’ 

Gj„ =c jl (0 x +c j 2 co y +c j 3 co z =[cjj c j2 C j3 ] ro y (2.b) 

_®z. 

where the subscript n denotes the nominal or design value. Combining all four gyros we obtain 



Fig. 2: The Gyro Configuration in the EOS-AQUA Satellite. 
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In the EOS-AQUA satellite the gyro configuration is as shown in Fig. 2 where, as mentioned before, X, Y, 
and Z are the axes of the body frame. The C matrix in this case is 


c = 



0 




(5) 


Note that the vector described by each row is of unit length like it should be (see Eq. 1). 

Due to misalignment, the sensitive axis of each gyro may actually point at a slightly different direction 
than the assumed one. This is illustrated in Fig. 1 where the components of this direction (which is still a 
unit vector) are cL, dj 2 , and dj 3 respectively. Following the steps that led to the development of the 
nominal gyro reading presented in Eq. (2.b), the actual gyro reading is found to be 

Gja = d ji©x +d j2 (Oy +d j3 co z = [dj, d j2 d j3 ] 


co. 


to. 


co. 


( 6 ) 


where the subscript a denotes an actual value. The difference, AGj , between the reading of the j — gyro 
and its assumed nominal reading is computable using Eqs. (2.b) and (6), as follows; 


AG” = G ja -G jn 




d j2 C j2 


dj3 C j3 j 


®x 

C0 y 

co. 


(7) 


where the superscript m denotes the fact that the error is due to misalignment. We denote by the 
differences d jt — c jj , dj 2 — Cj 2 , and dj 3 — Cj 3 as follows 


m jl= d jl- C jl (M m j2 =d j2- C j2 ( 8b ) m j3= d j3- C j3 ( 8 G 


Using Eqs. (8) we can write Eq. (7), as follows; 





m j2 


m 


j3 


(9) 


The m jj differences are shown in Fig. 1. Actually only two of the nijj of each gyro are independent. This 
results from the fact that the nominal as well as the actual directions of the gyros are given by vectors of 
unit length. For a reason that will be clear later, let us choose to present the third component of m j by the 
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first two; that is, we express m j3 , the misalignment along the sensitive axis, by m.j, and nij 2 . Since 
similarly to Eq. (1), it is also true that 

dji +dj2 + dj 3 = 1 0°) 

then using this relation and Eq. (8.c) we can write 

m j3 = ~ ^ji — dj2 — -y/l — c ji — c j2 GG 

For the case described by the fourth gyro (G 4 in fig. 3) the nominal direction of the gyro sensitive 
axis is along the Z axis; therefore, C 41 = C 42 = 0 , then from Eqs. (8.a) and (8.b) 

m 41 = d 41 (12.a) m 42 = d 42 (12.b) 

and from Eq. (11) 

m 4 3 =-\/l-d 41 -d 42 -1 (13) 

In the case where the misalignments are small, d 41 and d 42 are small too. Therefore we can expand the 
square root function of Eq. (1 1) in a Taylor series, as follows; 

>/l ~ d 41 - d 42 = 1 - yd 41 ~ 2"d 42 (14) 

(Note that the linear term of the series vanishes). Substituting of the last equation into Eq. (13) yields 

m 43 =-id 41 -±d 2 n (15) 

When d 41 and d 42 are indeed small, such as this case, then m 43 is negligible with respect to m 4] and 
m 42 . Then using Eqs. (12) we can write 


m 4i 

m 42 

m 43 


1 

0 

0 


Ml 


M2 


(16) 


It is the choice to express the component of ntj along the gyro sensitive axis (in this case m 43 ) that 
enables its elimination. 

For gyros whose sensitive axes are not aligned along one of the body axes the computation is 
more elaborate. Consider for example G 2 , the second gyro of the EOS-AQUA satellite. In order to define 
its misalignment in the body coordinates let us define a coordinate system in which the gyro sensitive axis 
is nominally aligned along one of its axes. Such a system (X”, Y”, Z”) is presented in Fig. 3, where the 

sensitive axis of the G 2 gyro is aligned along the system Y” axis. Following the preceding development 
for the G 4 gyro we conclude that 

m ' 2I = d 21 (16. a) m 23 =d’ 23 (17. b) 
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m 


22 


„ 2 ..2 

= -J-d -Id = 

2 U 21 2 U 23 


-i m ; 2 -^m ^ 2 


(17. c) 


(The “ sign denotes the fact that the values are expressed in the X”, Y”, Z” coordinate system.) Here too, 

•i 

the misalignment along the sensitive axis, m 22 , is normally negligible. In order to compute the 
misalignment error in the gyro reading we have to use Eq. (9) where the angular rate vector is transformed 
to the double prime coordinate system and the misalignment parameters are those given in Eqs. (16.a and 
b). As shown in Fig. 3 the transformation from the body to the double prime coordinates is performed by 


Z=Z’ 



Fig. 3: The Transformation from the G 2 Body to Gyro Coordinates. 

two rotations. The first rotation is by an angle a about the Z axis, and the second is by an angle — (3 
about the X’ axis. The resulting transformation matrix from the body to the G 2 coordinates is therefore 




ca 

sa 

0 



R-2 = 

-sa -cP 

ca-cP 

-sP 

(18) 



- sa • s(3 

ca-sp 

C p 


and 


to 2 = 

R> b 


(19) 

then following Eq. (9) 
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where 


It is easy to see that 


Define 


E 2 


AG;=[R^ b ] T d;=(oKd' 2t 

[«Q T =fe o d; 3 ] 


(20) 

( 21 ) 



ca 

-sa -sp 

_ 

R 2 b d' 2 , = 

sa 

ca-sp 

d 2l 

J " 


0 

cp 

_^23 _ 


( 22 ) 


ca 

-sa - sp 



a” 

sa 

ca-sp 

(23. a) 

d 2 = 

a 2 i 

A " 

0 

C P 



_d 2 3_ 


(23. b) 


then Eq. (20) can be written as 

AG™ = to^E 2 d 2 (24) 

To evaluate E 2 , we need to compute the angles a and (3 . For this we turn to Fig. 4 where these angles 


z 



Fig. 4: Definition of the Rotation Angles a and (3. 
are defined in the projection of the G 2 direction on the body axes. From this figure we conclude that 


a = tan' 


'21 


'22 


(25. a) 


P = COS 


-1 


,/c 7 " 


+ c 


22 


cos (y/ c 2 i + c 22 j 


(25. b) 


Using the EOS-AQUA satellite values (see Eqs. (4) and (5)) we obtain 
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a = tan 


-1 C 21 


= tan" 


- tan' 


P = cos' 1 V4 + c 2 22 = cos"' Uj + j = cos- 1 L/f 


= 35 . 26 ° 


(25. d) 


Due to the symmetry between the positioning of the G 2 and the G 3 gyros, it is easy to see that when 
considering the G 3 , the first rotation is about the new Z axis by the angle n — a and then about the X’ 
axis by the angle — (3 . Therefore 

- ca sa 0 

R 3 = -sa-cp -ca-cp -sp (26) 

- sa • sp - ca • sp cp 

and similarly to Eq. (19) for this transformation we obtain 

. . T> t> ^ /—I —I \ 


then following Eq. (20) 
where 

It is easy to see that 


®a =R 3 ca b 


AG™ =[R 3 b <o b ] T d; =<R- b d; 


Define 


[d;.f=[d 

o 4 ] 

-ca 

-sa - sp 

R b d 3 * = sa 

-ca-sp 

0 

cp 

- ca - sa • sp 


sa - ca ■ sp 

0 cp 

(31. a) 


(3 1 b) 


then Eq. (28) can be written as 


AG™ =co 1 E 3 d 3 


For the G 3 gyro we have only one transformation, which brings the body X-axis into coincidence 
with the Gj gyro sensitive axis (see Fig. 2). It is about the Y-axis by an angle which we denote by y . For 
this gyro we have then 

cy 0 - sy 

R b = 0 1 0 (33) 

sy 0 cy 


and the angular rate in coordinate system 1 is then 
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(O, =R 1 b e> b 

(34) 

and 


AG”=[R> b ] T d;=toXd;. 

(35) 

where 


[d;,] T =[0 d u d 23 J 

(36) 


We denote the misalignment parameters of this gyro by a single prime because it takes only one rotation (to 
a single prime system) to align the coordinate axis with the sensitive axis of the G, gyro. Note that here 

the misalignments that are not negligible are d 12 and d 13 . It is easy to see that 


Define 


E, 


then Eq. (35) can be written as 


RUn = 


0 sy 

1 0 
0 cy 


*12 


'0 

sy 



r n 

1 

0 

(3 8. a) 

d,= 

dn 

0 

cy 



L d 13 J 


(38.b) 


AG” = ct>bEjd, 


(37) 


(39) 


From fig. 2 it is easy to see that the rotation angle, y , is computable, as follows; 


y = sin 


-i 


•'ll 


= sin 


1 1=35.26° 


(40) 


Similarly to the computations carried out for the misalignment errors for gyros 1, 2, and 3, we can 
write for gyro 4 

AG;=<E 4 d 4 (41) 

where, based on Eq. (16), 



1 

0 





e 4 

= 0 
0 

1 

0 

(42. a) 

d 4 = 

^41 

41 (42. b) 

_ 42 _ 


[ag” 

r-i 

AG” 

AG” 

AG” 

ag:] 

(43. a) 


Q m 


0 

0 

0 


CD y co 2 0 0 0 

0 o co x 0) y G) z 

0 0 0 0 0 

0 0 0 0 0 


0 0 

0 0 

®x ®y 

0 0 


0 0 0 0 

0 0 0 0 

© z 0 0 0 

0 co x CD y co z 


(43. b) 
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E = 


0 


0 E, 


0 

0 


0 

0 


0 E, 


0 

0 

0 

E, 


(43. c) d T = [d* d 2 dl d 4 r ] (43. d) 

U | 

0 0 

Then Eqs. (39), (23), (32) and (41) can be unified into the following single equation 


AG m = Q m Ed 


(44) 


Scale Factor Error Model 


As mentioned, another error source that causes the difference between the correct value of the 
rates and their measurements are the scale factor errors. The error model for the scale factor error is simply 


AG k 


m G1 k, 

®G2 k 2 

G>G3 k 3 

®G4 k 4 


(45) 


where the subscript k denotes the fact that this error is caused by gyro scale factor error, C0 Gj , i = 1 - 4 is 

the angular velocity measured by gyro number i, and k ; is the scale factor error of that gyro. The actual 

components C0 Gi are obtained by transforming the angular velocity expressed in body coordinates to the 

actual misaligned gyro sensitive axes using the matrix D; however, because D is unknown to us we use 
instead the matrix C that transforms the body rate to the nominal gyro axes, and is close enough to D. Thus, 


®G1 


_ 



co x 

®G2 

= c 

tOy 

®G3 


y 



_®z_ 

_ ra G4_ 




(46) 


where in the case of EOS-AQUA, C is as given in Eq. (5). Equation (45) can be written as follows 


Define 


1 

s 

o 

0 

0 

1 

O 

V 

0 

0) G 2 

0 

0 

k 2 

0 

0 


0 

k 3 

1" 

o 

0 

0 

®G4_ 

_ k 4. 


Q k 


®G1 

0 

0 

0 


0 0 0 

cd G2 0 0 

0 co G3 0 

0 0 ro G4 


(47) 


(48. a) 
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and 

k T = [kj k 2 k 3 k 4 ] (48.b) 

then Eq. (47) can be written as 

AG k =Q k k (48. c) 

Bias Model 

The bias error model is quite simple and is given by 

V 

b 2 

b = 2 (49) 

_ b 4. 

where b, is the bias of gyro number i. 

The Augmented Gyro Error Model 

The total gyro error is the sum of all the errors discussed before; namely, bias, scale factor and 
misalignment errors; that is 

AG = AG m + AG k + b (50.a) 

or using Eqs. (44) and (48. c) 

AG = Q m Ed + Q k k + b (50.b) 

The last equation can be written in the following form 

Id' 

G a -Cco r = [q“E Q k I 4 ] k (50. c) 

b 


where 0) r is the reference angular velocity vector. It is the angular velocity, which the SC experiences in 
reality. As mentioned before, G n is the nominal angular velocity measured by the four gyros. The vector of 
the left hand side of the last equation as well as the matrix on the right hand side are functions of the body 
angular rate, (O r . We denote them as follows 


y(<*>r) = G a - Ct »r 


(51. a) 


also let 


then Eq. (50.c) can be written as 


H(to r ) = [q™E Q k 



y(oj r ) = H(<o r )x 



(5l.b) 
(51c) 
(5 Id) 
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CALIBRATION-PARAMETERS ESTIMATION FOR KNOWN RATE 


Our goal now is to estimate X , and for that we need to know the angular rate, which the gyros are 
set to measure. We distinguish between two major cases. One case is that where tO r , the reference SC 
angular velocity, is known, and the other case is that where the rate is not known and has to be evaluated 
simultaneously with the estimate of X . In the first case we also distinguish between the deterministic and 
stochastic cases. All these cases are discussed next. 

Deterministic Case 

When the SC rotates at a certain angular rate and a one-time measurement of the four gyro readings is 
taken at that time, which we denote by t k , we obtain one matrix equation, as follows: 

y(®r,k) = H (<°r,k) x (52) 

where CO r k denotes the angular rate at time t k . This yields four equations for the 16 unknowns of X . If the 

rate does not change, then more measurements do not change the equations. A change in the angular rate of 
the SC is needed to generate more equations. It should be noted that even if we have 16 equations, it does 
not mean that they are all independent and that we can solve for X . We have to design the profile of 

co r and the times when measurements are taken in such a way that we will be able to find 16 independent 
equations. Let us denote the 16 independent equations by one matrix equation, as follows: 

y = Hx (53) 

Because we have 16 independent equations, H has an inverse; therefore, we can solve for X using 

x = H“'y (54) 


Stochastic Case 

In this case we assume that the measurements are contaminated by noise, which is the most likely 
case. Therefore, the matrix equation that describes this case at time t k is 

y( W r,k) = H (®r,k) X + V k ( 55 ) 

Even if we find 1 6 independent equations from measurements done at different time points we still want to 
use all available measurements and obtain x as a least squares estimate. We have 


’y(®,.)’ 


H(o rl )~ 


’Vj" 

yK.z) 


H(e> r , 2 ) 

x + 

V 2 

yO r , 3 ) 


H(G> r , 3 ) 

V 3 

jyK.J 


H(o r>4 ) 


_ V 4 _ 


Let 



(57 .b) 


y(°>r,l) 


H(® rl ) 

y(<°r, 2 ) 

II 

DC 

'cd' 

hk >2 ) 

yO r , 3 ) 

H(o> r3 ) 

_y(«r,4>_ 


H(© r>4 ) 


then X , the least squares fit to X , is as follows (ref. 5) 

x = (H t H)“'H t Y (58) 

The profile of co r has to be chosen in a careful way as to enhance the observability of X . 


CALIBRATION-PARAMETERS ESTIMATION FOR UNKNOWN RATE 

In this case we have to find the angular rate vector while estimating the calibration parameters. The 
information that we have is attitude information and gyro measurements. We need the attitude information 
in order to estimate the angular rate, and we need the gyro measurements, as well as the estimated angular 
rate, for the calibration process. The attitude information can be supplied in various ways; namely, we may 
have it in the form of raw vector measurements or we can have it in an already processed form as attitude 
quaternion for example. The angular rate behaves according to the following SC angular dynamics equation 

d) = I -1 [(1(0 + h)x]to + 1 -1 (T - h) (59. a) 


where I is the SC inertia tensor, [(Ico + h)x] is the cross product matrix of the vector (Ito + h) , h is the 

angular momentum of the momentum wheels, and T is the external torque operating on the SC. Because x 
is a constant vector, it obeys the following differential equation 

X = 0 (59.b) 

We are tempted to combine the last two equations into one dynamics equation, as follows; 


O) 


~r‘[(Ito + h)x] 

o' 

0)1 

-f 

I _i (T - h) 

X 


0 

0 

! 

X 


0 


This dynamics model calls for the use of a Kalman Filter (KF). In fact the most appropriate filter is the 
Pseudo-Linear Kalman (PSELIKA, ref. 6) filter. To find the suitable measurement equation, we turn to Eqs. 
(50. c) and (51.b) from which it is obvious that 


G, =[C 



(60) 


In order to apply the filter algorithm we need to add some process noise to Eq. (59. c) and some 
measurement noise to Eq. (60). However, it is easy to see though that the last system is unobservable. 
Theoretically, if the actual noise values are small, and if we know the initial angular velocity, then we can 
compute the angular velocity separately. Once we know the angular rate and we command the SC to 
execute a suitable angular rate profile, we should be able to compute the calibration parameters. In reality 
though the computation of the angular rate is not accurate enough because it is done in an open loop 
manner. Therefore, we have to add attitude measurements in order to check the divergence of the computed 
angular rate. In this case we can indeed combine all the dynamics equations into one augmented matrix 
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equation and estimate the augmented state vector. This is so because it is possible to find an angular rate 
profile that will render the system observable. 

As mentioned, attitude can be given in several ways; namely, it can be given in a raw form as 
vector measurements or in processed attitude parameters like a quaternion or direction cosine matrix 
(DCM). Let us consider two cases, one where attitude is represented by a processed quaternion, and the 
other case when we have raw vector measurements. The case where attitude is given in the form of a DCM 
can be inferred from the development presented in Ref. 7 and the way we handle quaternion representation 
of attitude. 

Estimation When Attitude is Presented by the Attitude Quaternion 

Let us assume first that the attitude is given in a form of a quaternion (ref. 8). In this case the filter 
dynamics is as follows (ref. 7) 


o') 


I -1 [(Ito + h)x] 

0 

o’ 

CO 


’r'(T-h)’ 

X 

= 

0 

0 

0 

X 

-1- 

0 



IQ 

0 

0 

_q_ 


0 


where 


q 4 

-q 3 

q 2 

q 3 

q 4 

-q, 

-q 2 

qi 

q 4 

-q. 

-q 2 

-q 3 


and the corresponding measurement equation is 

fto 



(61) 


(62) 


(63) 


The matrix I 4x4 is a fourth order identity matrix. The combined measurement equation consists of Eqs. 
(60) and (63); that is 


G a 


' c 

H(«) 

0 4 x 4 

_q,n_ 


_0 4 x3 

^4x16 

Qx 4 _ 


(64) 


Estimation When Attitude is Given by Vector Observations 

Normally, in space missions attitude is determined from vector observations. These observations can be 
used directly to check the divergence of the angular velocity estimates (ref .7). This is shown next. Suppose 

that we have N vector measurements at a certain time point. Let r ; denote some abstract i* vector as 

expressed in the reference coordinate system, and let denote the same vector when expressed in the body 
coordinates. From the laws of dynamics it is known that 

Dr ; = b; +coxb j (65) 
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where r ; is the time derivative of r { as seen by an observer in the reference coordinates, D is the matrix 
that transforms vectors from the reference to body coordinates, and b, is the time derivative of b, as seen 

by an observer in body coordinates. The vector b ; is a measured vector and b s is its time derivative. We 
can write Eq. (65), as follows; 

b, =[b i x]c» + Dr i (66) 

Note that r ; is computable since 1* is usually known because, generally, the vector is a direction to a 
certain known planet whose location is given in an Almanac or, like with magnetometer measurements, the 
vector can be computed using a model. (It should be noted that quite often the rate of change of r is so 
small that r, is negligible). Define 


V 


[b,x] 


Drj 


(67.a) B = 


(67. b) and U = 


>n. 


_[b N x] 


_ d *n_ 


(67 .c) 


then we can augment all the N equations of Eq. (66) into one matrix equation, as follows; 


P = Boo +u 


( 68 ) 


therefore, instead of Eq. (61), we obtain in this case of vector measurement the augmented equation 


to 


'r’Kio+^x] 

0 

o' 

to 


’r'(T-h)' 

X 

= 

0 

0 

0 

X 

+ 

0 

P_ 


B 

0 

0 

_P_ 


u 


and the corresponding measurement equation is 

Pm = |p 3x3 ^3x16 


whereas the augmented measurement equation is 



'G. 

Pm 


C U(to) 0 


0 


3x3 


4x4 


0 


3x16 


I 


3x3 


(69) 


(70) 


(71) 


COMPENSATION 

To complete the calibration process we need to perfonn its second stage; namely, compensation 
where we eliminate the estimated errors from the gyro readings. From Eq. (50.c) we obtain 
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Coo r = G 


(72. a) 


a -[Q m E Q k 



d 

k 

b 


As mentioned before, G a is a vector of the gyro readings and (D r is the correct angular velocity vector. 
However, we do not have CO r which is what we are trying to measure; therefore, to compute Q m and Q k , 


which have to be computed using (0 r , we use the measured uncompensated angular rate vector. This 


vector is derived from the uncompensated gyro measurements which we called actual and denoted by a. 
Also, we do not have the actual values of d, k or b, but rather their estimate; therefore, using the values on 
hand Eq. (72. a) becomes 


C<b r 


g,-[o;*e n; 



<72.b) 


where a denotes the actual values and A denotes estimated vectors. To obtain the compensated 
measurements of the angular rate vector define the A matrix, as follows: 


It is easy to verify that 


A = 


1 0 
0 1 
0 0 


0 

0 

1 

2 . 


ac t c = i 


therefore pre multiplying Eq. (72. b) by AC r yields 


<a r = AC T G a -AC T [Q a m E Q k 



(73) 


(74) 


(75) 


SIMULATION RESULTS 

In lieu of actual SC data, a simulator was developed to produce the gyro, reaction wheel, and star 
tracker data. Much care was devoted to the simulation since any dynamics errors would be perceived by 
the KF as a state error. First, the maneuver strategy was developed. An inertial period before any 
maneuver would facilitate the estimation of the gyro bias. The scale factors of each gyro could be 
estimated by a maneuver about that gyro axis. The misalignments could be estimated by the same scale 
factor maneuvers. An additional two maneuvers about the SC X and Y body axes, respectively, were added 
to assist in the alignment estimation of gyro 4 which senses rate about the body Z axis. Second, the 
maneuvers were modeled as a ramp-coast-ramp where the linear ramp time was 5 seconds. The rate profile 
was then obtained and CO was derived by simple subtraction. Third, the dynamics Eq. (59.a) was re-written 
in terms of hand the ordinary differential equation (ODE) was solved using the to and CO generated 
above. Fourth, using the newly generated system momentum profile h , the rate, CO , was determined from 
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Fig. 5: Gyro Bias Estimate (bold) versus Truth (thin) 



Time (sec) 

Fig. 6: Gyro Scale Factor Estimate (bold) versus Truth (thin) 
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Fig. 7: Gyro Misalignment Estimate (bold) versus Truth (thin) 


the ODE in Eq. (59.a). Fifth, using this 0) , the quaternion was estimated using the q portion of Eq. (61). 
Lastly, the body rate was resolved in the AQUA gyro frame by use of Eq. (50.c). 

The KF was then executed on the simulated data. The biases were estimated using the initial 
inertial period. Also, to facilitate the bias estimation, the scale factor and misalignment parameter 
estimation was terminated by zeroing out their respective rows and columns in the state covariances and 
process noise. The scale factors for each gyro were then estimated using the respective maneuver about 
that axis and by zeroing out the influences of the other scale factors, the misalignments, and biases as 
described with the bias estimation. Lastly, the gyro misalignment slews, which were a repeat of the scale 
factor slews with the addition of X and Y-axis maneuvers, were executed with the biases and scale factors 
zeroed out. All states were estimated with less than a 1% deviation from truth. The KF bias estimate can 
be seen in Fig. 5, followed by the scale factor estimate in Fig. 6, and lastly the misalignment estimate in 
Fig. 7. 

CONCLUSIONS 

In this paper we presented a new method of gyro calibration. Normally, we have to calibrate a 
cluster of three gyros whose sensitive axes are along the body axes. Here, the rate is read by four gyros only 
one of which is aligned along the body coordinate axes. Therefore, a new algorithm was devised for 
calibrating a quadruplet rather than the customary triad gyro set. In particular, a new model had to be 
developed for the gyro misalignment errors. Normally, the gyro outputs are used to supply data for a 
differential equation, which is solved in order to compute attitude. According to the new method the gyro 
outputs are also used as measurements, which are fed into a Kalman filter that estimates the gyro 
misalignment, scale factors, and biases. The new calibration algorithm was developed in particular for the 
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calibration of the EOS-AQUA satellite gyros. The effectiveness of the new algorithm was demonstrated 
through simulations with error of each estimated parameter being less than 1%. 
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